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ABSTRACT 

Fermi-LAT analyses show that the 7-ray photon spectral indices T^ of a large sample of blazars 
correlate with the uJG peak synchrotron frequency Vg according to the relation T^ = d—k log Vg. The 
same function, with different constants d and fc, also describes the relationship between T-y and peak 
Compton frequency vq- This behavior is derived analytically using an equipartition blazar model 
with a log-parabola description of the electron energy distribution (EED). In the Thomson regime, 
k = kEC = 36/4 for external Compton processes and k = kssc = 96/16 for synchrotron self-Compton 
(SSC) processes, where 6 is the log-parabola width parameter of the EED. The BL Lac object Mrk 
501 is fit with a synchrotron/SSC model given by the log-parabola EED, and is best fit away from 
equipartition. Corrections are made to the spectral-index diagrams for a low-energy power-law EED 
and departures from equipartition, as constrained by absolute jet power. Analytic expressions are 
compared with numerical values derived from self-Compton and external Compton scattered 7-ray 
spectra from Ly a broad-line region and IR target photons. The T.^ vs. Vg behavior in the model 
depends strongly on 6, with progressively and predictably weaker dependences on 7-ray detection 
range, variability time, and isotropic 7-ray luminosity. Implications for blazar unification and blazars 
as ultra-high energy cosmic-ray sources are presented. Arguments by Ghisellini et al. (2014) that the 
jet power exceeds the accretion luminosity depend on the doubtful assumption that we are viewing at 
the Doppler angle. 

Subject headings: gamma rays: galaxies—radiation mechanisms: nonthermal—galaxies: jets— 
galaxies: BL Lacertae objects: general—galaxies: quasars: general—acceleration 
of particles 


1. INTRODUCTION 

Searches for an ordering principle in blazar science 
have met with limited success. One of the most de¬ 
bated is the blazar sequence, in which blazar data seem 
to show an inverse correlation between apparent isotropic 
synchrotron luminosity Lgyn and peak synchrotron fre¬ 
quency Vg of the blazar z/JC spectral energy distri bution 
(SED) (jFbssati et al.lll99^ ISambriina et al.lll996[l . This 
behavior, which is mirrored in the 7-ray r egime, has been 
interpreted in terms of cooling processes (iGhisellini et al.l 
119981 : iBbttcher fc Dermeill20n^ lFinkel[^l,'lll . The valid¬ 
ity of the blazar sequence h as, however, been criticized 
(|Giommi et al.l 1201^ l2013ll as possibly resulting from 
spurious correlations introduced by combining samples 
from radio and X-ray blazar surveys, problems from red- 
shift incompleteness, and confusing lineless BL Lac ob¬ 
jects that lack accretion disk with those where the BLR 
radiation is overwhelmed by beam ed emission. Con¬ 
trary to the simple blazar sequence, iMever et al.l ()2nil[ l 
present evidence for the existence of two separate tracks 
in the Lgyn vs. Vg plane, including radio galaxies in the 
blazar-sequence plot. 

A second strong correlation is the blazar divide. From 
the first three months of Fermi Large Area Telescope 
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(LAT) blazar data, IGhisellini et al.l (|2009ll argued that 
hard (F.^ < 2) 7-ray spectrum blazars are associated 
with sources radiating isotropic 7-ray luminosities Lj < 
5 X 10®® erg s“®, while soft (F.^ > 2) 7-ray blazars are 
more likely to be at larger values of From the Sec- 
ond Fermi LAT AGN (2LAC) data (lAckermann et ahl 
Enni), a broad divide is evident in t he direct data at 
= 10®® erg s“® (Figs. 37 and 38 in lAckermann et al.l 
I2OIIII , though with no other apparent dependence of F.^ 
on Lj in the ranges 10®® ^ L.y < 10®® erg s“® and 
10®® ^ L7 ^ 10®® erg s“®. In terms of a beaming- 
corrected Eddington ratio ^Edd for a black hole with mass 
^ 10® Mq, this could imply a transition from an ineffi¬ 
ciently radiating ADAF-typ e flow at fudd 0-0 1 to a 
thick disk when ^Edd ^0.01 (jChisellini et al.l[2009ll . The 
7-ray Compton dominance Ac, which is essentially the 
ratio of the bolometric 7-ray and sy nchrotron luminosi¬ 
ties, also st rongly correlates with Vg (jFossati et al.lll99^ 
iFi^ [201311 . 

Definitive interpretations of blazar sequence and blazar 
divide data are hampered by redshift incompleteness. BL 
Lac objects without redshift information may themselves 
constitute separate populations in the Lgyn vs. Vg or F.^ 
vs. L7 planes, though large efforts have been made to 
provide complete, or at least redshift-constrained sam - 
ples of blazar data (IShaw et al.ll2013l : lAiello et al.ll2014ll . 
The Ac vs. Vg distributions of 2LAC blaz ars with an d 
without redshift do not significantly differ (|Fin SiMl). 

A third robust correlation in blazar physics relates 7- 
ray spectral index F-j, with peak synchrotron (vg) or peak 
Compton 7-ray (uc) frequencies (in this study, we as- 
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sume that the blazar SEDs are made by leptonic pro¬ 
cesses only).® These spectral-index diagrams for FSRQ 
and BL Lac blazars have been re ported in the First LAT 
AGN Catalog (ILAC, Fig. 13 in lAbdo et al.l|2010all. the 
Fermi Bright Blazar SED pape r (Fig. 29 in|AbciQ_gL^ 
I 2010 bli . the 2LAC (Fig. 17 in lAcker mann et al.l l 2 nilli . 
and the 3LAC (Fig. 10 in lAckermann et al.ll2015li . The 
distributions of spectral indices of the entire BL Lac 
and FSRQ blazar samples follow a pattern, with large 
scatter, described by the relation T..^ = d — A: log 1 ^ 14 , 
where i/g = Hz. For the entire sample of FS- 

RQs and BL Lac objects, t he va lue k = 0.18 ± 0.03 is 
found in lAckermann et al.l (|2015h . A similar function, 
with different values of d and k, apply to the L.^ vs. 
i^c data. The spectral-index distribution of BL Lac ob¬ 
jects with unknown redshift is generally consistent with 
the distribution of BL Lac ob jects with known redshift 
(|Ackermann et ^l201lL l2015ll . 

In this pape r, we use an equip a rtition blazar model- 
ing approach (|Cerruti et al.l l2013l : iDermer et al.l l2014a^ 
bf assuming a log-parabole electron energy distribution 
(FED) to explain the blazar spectral-index diagrams. In 
Section 2 we derive analytic Thomson-regime expressions 
for the relationship between L.^ and Vg, depending on 
whether the 7 rays are made through external Compton 
(EC) or synchrotron self-Compton (SSC) processes. Be¬ 
cause of the equipartition relations, the expressions de¬ 
pend on b, j/g, variability time tyar, bolometric isotropic 
synchrotron luminosity Lsyn, an equipartition parame¬ 
ter and a radiation parameter Cg. The simpler T.^ vs. 
VC expressions are also obtained. The derived analytic 
relations, confirmed by numerical modeling, are shown 
in Section 3 to be in general accord with the blazar 
spectral-index diagram data, whether external radiation 
helds in the jet environment are present or absent. The 
effects of a log-parabola FED with a low-energy power- 
law component are also considered. 

In Section 4, application of the equipartition model 
to the BL Lac object Mrk 501 is demonstrated, and 
effects of departures from equipartition are evaluated. 
Trends in spectral-index behavior with other observables 
constrained can be tested with correlated Fermi-LAT 
and multiwavelength data, and how this work relates to 
the blazar sequence and blazar divide, and blazars as 
UHECR sources, are discussed in Section 5. The work 
is summarized in Section 6 . 

Appendix A gives a Thomson-regime derivation of the 
SSC spectrum with a log-parabola electron distribution, 
and Appendix B gives a jet power-analysis. There 
we show that the assumption that blazars are typi- 
call y observed at the Dop pler beaming angle may have 
led IChisellini et all (| 201 ^ to overestimate the absolute 
jet power. Indeed, out-of-equipartition models are ulti¬ 
mately constrained by demands for power. 

2. EQUIPARTITION BLAZAR MODELING WITH 
LOG-PARABOLA ELECTRON ENERGY DISTRIBUTION 

A standard blazar-jet model, treate d in innumer¬ 
able blazar spectral modeling papers (see iBdttcher et al.l 
I 2 OI 2 I for review), starts with magnetized plasma that is 


ejected at relativistic speeds along the poles of a rotat¬ 
ing black hole.® The jet plasma, which entrains thermal 
and nonthermal particles in a hypothetical tangled and 
randomly oriented magnetic field, is a source of escaping 
photons, and potentially also of escaping cosmic rays and 
neutrinos. The jet power is extracted from the mass en¬ 
ergy of accreting matter and/or the rotational energy of 
the black hole itself. The collimated relativistic plasma 
outflow, an exhaust byproduct of the energy generated by 
the black-hole engine, is usually attributed to processes 
taking place in the magnetosphere of the rotating black 
hole. The polarized broad-band synchrotron radiation 
emitted by an energetic FED (which could also contain 
positrons) is boosted by the Doppler effect along the jet 
axis, so that rapidly variable jet synchrotron radiation 
can be detected by Earth-based observatories from large 
redshift {z » I) sources. 

The jet electrons also Compton scatter ambient pho¬ 
tons to 7 -ray energies. Besides the accompanying 
SSC emission from t arget synchrotron photo ns (e.g., 
iMaraschi et al.l 119921 : iBloom fc Marsch^ 119961) . EC 7 
rays are made when the nonthermal jet electrons scat¬ 
ter photons from external radiation fields. Depend¬ 
ing on jet Doppler factor (5 d and BLR cloud param¬ 
eters, the direct accretion-disk radiation field domi¬ 
nates the external radiation field of a powerful FSRQ 
at <C 10® Schwa rzschild radii, BLR fields are strongest 
within ^ 0.3 pc (IDermer et ai.lll99l : ISikora et al.l 11994 
IDermer fc Schlickeised 120021 ). while at the pc scale and 
beyond, infrared radiation from a surrounding IR- 
emitting dust torus would have t he largest energy density 
of all ambient rad i ation fields (iB lazeiow ski et al.l 120001 : 
ISikora et al.l 120091: IChisellini fc Tavecchiol 120091 ) in the 
inner blazar-jet environment. 

The form of the nonthermal FED is often treated by 
either assuming a nonthermal injection spectrum of lep¬ 
tons that evolves in response to adiabatic and radiative 
losses, or by assuming a form for the average steady-state 
FED in the radiating jet plasma. Adopting the latter 
approach, we assume that the 3 parameter log-parabola 
function 

—61og(^7—) ,, 

^pk 

provides an approximate description of the nonthermal 
lepton spectrum. Here y = 7V7pfcj 7pfc peak, or 

principal, Lorentz factor of the fluid-frame EED, eq. O- 
The value of K' can be related to either th e total particle 
number or total comoving particle energy (IDermer et al.l 
1201 4al) : in the latter case, K' = In 10/6, 

where is the nonthermal electron energy of the blob. 

The continuously curving EED given by a log-parabola 
function derives from stochastic acceler ation processes 


with radiation and escai 

oe (see, e.e.. IMassaro et al.l 

200l iBecker et al.l IMiS 

Tramacere et al.l l2007l 120111: 

Stawarz & Petrosianll2008D 

. With this form of the EED, 


GeV breaks in FSRQs and blazars with 1 /®/" < 10®^ Hz 


are shown to arise from the onset of Klein-Nishina ef- 
fects when scattering BLR photons (iCerruti et al.ll20f^ 


y^A^'(7') = {I'^kK^i'pk 


® A correlation of the log-parabol a width pa rameter b and Us 
is apparent in SED modeling studies l|Chenll2014l') , but is based on 
only 5 or 6 high-synchrotron peaked blazars. 


® The shock-in-jet model of IMarscher Geari II1985D provides an 
alternate approach that could apply to the 10^^ Hz radio regime 
that often remains unfit in the standard model described here. 
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TABLE 1 

Dependences of 5d, B ', and 7]'^^.“ 



Coef. 

A48 

1/14 

u 

Cs 

Ce 

/o 

/i 

/2 


17.5 

3/16 

1/8 

-1/8 

-7/16 

1/4 

-7/16 

-1/4 

-1/8 

B'{G) 

5.0 

-1/16 

- 3/8 

-5/8 

13/16 

-3/4 

13/16 

3/4 

3/8 

'^pk 

523 

-1/16 

5/8 

3/8 

-3/16 

1/4 

-3/16 

-1/4 

-5/8 

e‘> 

1.4 

5/16 

- 1/8 

1/8 

-1/16 

-1/4 

-11/16 

1/4 

1/8 

T c 

4 

5/8 

- 1/4 

1/4 

-1/8 

-1/2 

-1/8 

1/2 

1/4 


“ So, e.g., 5 d ^ 17.5L®/^®(l^l4//2^4)l/«(/oC.)■^^^®(Ce//l)l/^ 
etc. 

^ S — Emax (1020 eV)/Z 

^ Absolute power in magnetic field, units of 10^“^ erg s“^ 


lAckermann et aT I I2OIOI: iTavecchio fc Ghisellinil 120081) . 
and to give ( Dermer et al.ll2014al) reasonable fits to four 
epochs of quasi-s imultaneous multiwave length observa¬ 
tions of 3C 279 (|Havashida et al.ll2012i) . As we show 
below, this approach also gives good fits to the SED of 
Mrk 501, though the best fits are achieved with an elec¬ 
tron distribution out of equipartition with the magnetic 
field. 

The comoving synchrotron luminosity 


n/2 roo 

L'syn = CTT-^ ( 2 ) 

implies, using Eq. o and a (5-function approximation 
for the synchrotron photon with ave rage dimensionless 
energ y Cgyn = /Bcr)"/'"^ (jPermer Sz Menoiil 

I2OO9I ). the received synchrotron luminosity spectrum 


eLsynie) = = t;( —, (3) 

^pk 


where x = ^/e / tpk, v = hLsyn, a nd = 2 • 

In 10/6 (jPermer et al.l I^Oldal) . Thus the ef¬ 
fective log-parabola width parameter bgy for the syn¬ 
chrotron spectrum is given by bgy = 6/4 in the 5- 
function approximation, and bgy = 6/5 when using 
the full Thomson cross section (|Massaro et al.l 120061 : 
iPaggi etld] I2009D . The peak synchrotron frequency 
Cpfe = {2>l2)5Yy{B’/Bcr)lpk- The slope of the eLsynie) 
spectrum is 


=l|i-i,iog(tA,D]. (4) 

(|Massaro et al.l l2004aD . Because the nonthermal elec¬ 
tron energy-loss rate from synchrotron processes scales 
quadratically with electron Lorentz factor 7 / the syn¬ 
chrotron spectrum from a log-parabola distribution of 
electrons has a synchrotron peak energy Cg = 

hvslme,(? that is shifted to higher values than tpk- Eq. 
(HI shows that = lO^/^Cpfe (|Massaro et al.l[2?)0^ . 

Table[T]shows the various dependencies of blob proper¬ 
ties on the observables L 48 , Cs (or 2 ^ 14 ), tyar and 6 , and on 
the equipartition factor Ce and radiative factor (/g. The 
factor Ce is the ratio of nonthermal electron energy den¬ 
sity itg to magnetic-field energy density u'^, = B'^/Stt, 


and Cs is the ratio the jet-frame synchrotron photon en¬ 
ergy density and u'g,. In the blob scenario, the geome¬ 
try factor /o = 1/3. The 6 -dependent factors are /i = 
10 -i/ 4 fc, fy = 10 V^ and /s = (2 • 10/6)-4 

(|Permer et al.ll2M4al) . so Cg = f 2 ipk- 
The similarity of the underlyin g physics of the syn¬ 
chrot ron and Compton processes (|Bhimenthal fc Gouldl 
means that an expression like eq. 0 holds for 
Compton scattering in the Thomson regime, except now 
Cpk is replaced by a correspon ding peak photon energy 
for EC and SSC processes ('e.g.. lPaggi et al.ir2009D . In the 
EC case, epk,EC = ( 4 / 3 ) 6 Deo 7 /fc, assuming an isotropic 
monochomatic external radiation field with energy eg 
and energy density u n- From the equipartition relations 
(jPermer et al.l[^I4al) shown in Table I for the photon 
spectral index P.^, = 2 — 0 :^, we find that the photon index 
for EC processes is given by 


^C = w + ^og 


1 1/2. .1/4 

\^Lya^4 Se4^48 , 


- ^ log 1414 . (5) 


Here BceV is the effective detection energy in GeV, and 
Co = 2 X 10 ~^eLya for Ly cr/BLR scattering. A nominal 
value of EoeV = 1 is chosen because the Fermi-LAT is 
most sensitive at ~ 1 GeV (Fig. 18 in lAbdo et al.ll20104 
for a r.y = 2.2 source spectrum). The dependence of T..^ 
on EoeV can be studied by analyzing Fermi-LAT data in 
discrete energy ranges. 

Scattering the dusty torus emission, with IR pho¬ 
ton energies corresponding to e^ya 0 . 02 , implies a 
Thomson spectrum softer by AT-^ = 0.856, because 7 - 
ray photons at a given observing energy are produced 
in the softer part of the Gompton-scattered spectrum 
when the target photons have lower energies. If equipar¬ 
tition is instead made to total particle energy density 
according to the factor Qg = u^g^jug,, then Q = 
Ceq/{^ + Vbi) and r]u = Uy^ryonsl'^e Is the baryon loading, 
and le the i nternal energy density in protons and 

ions (iDermer et al.ll20i4H and Appendix B). 

The specific spectral synchrotron luminosity, from 
eq. (H) in the ^-func tion approximation and results of 
IDermer et al.l (l2014aD . is given by 


eL,,„(e, H) = /gV, ^ caxl^ , (6) 


where x = and tpk = 4(5Di?'7p|/3i?cr- The spe¬ 

cific spectral 7 -ray luminosity in the Thomson regime for 
a jet traveling through an external isotropic, monochro¬ 
matic radiation field with frequency m1eg/h and energy 
density ug, in units of rrieC^ cm“^, using a ^-function 
approximation for Thomson scattering, is 

eLgcie, n) ^ hNg ^ ca^ug ^ (7) 


where v = ^ Jejapk^EC and Cpk^EC = ( 4 / 3 ) 6 g) 7 ((|eo. The 
technique of iGeoreanopoulos et al.l (l2001l) is used to de¬ 
rive this expression. The ratio of the spectral syn¬ 
chrotron and Thomson luminosities at their respective 
peak frequencies is S'^ug/u'g,. 

For the SSG process, the combined effects of the widths 
of both the EED and the target synchrotron photon 
spectrum will broaden the Gompton-scattered photon 
spectrum such that its effective width in the Thomson 
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regi me is obtained by replacing b by bssc = b/2 in eq. 
(gl) (jPaggi et a ll 1200911 and replacing Cpfc by epk,ssc = 
2dDiB'/Bcr)%\, giving 

T.55C ( 6.5 X lO^EGevfo^\s^^Lli^ \ 

i 7 32 ^ 4 ® ^ ) 

- ^ log 1^14 , (8) 

from Table 1. Note that the lyL^, peak SSC frequency is 
a factor 10^/^ larger than epk,ssc- The SSC expression 
is justified by a more detailed derivation in Appendix 
A. The uncertainty AT^ in the spectral index related 
to geometrical uncertainties can be estimated by letting 
/o range from unity for a blast-wave shell geometry to 
/o = 1/3 for a comoving spherical-blob geometry. From 
eqs. (ED and ED, one can see that this translates into 
an uncertainty = 0.306 for EC processes and an 

uncertainty AF®'^'" = 0.046 for SSC processes. 

We also derive the Thomson-regime expressions 

F^^’T' = 2 + ^ log(2.4F;Gev) - I log 1^23 (9) 

for F.^ vs. VC in EC processes, and 

rSSC .7 = 2 + ^ \og{2AEGev) - 2 log ^23 (10) 

for r.y vs. VC in SSC processes. Here 1^23 = i^c/10^^ 
Hz is the peak frequency of the Compton component of 
the vLi, SED. Note that eq. ED is independent of the 
target photon energy, because the expression assumes 
that the EED and Doppler factor are adjusted to produce 
a Compton-scattered 7 -ray spectrum that peaks at vq. 

3. MODELING THE BLAZAR SPECTRAL-INDEX 
DIAGRAM 

Fig. [T] shows measured values of Fermi-LAT sp ectral in¬ 
dex F.^ from the 2LAC (jAckermann et aHl2011ll derived 
from a single power-law fit to the complete data set in 
the 0.1 - 100 GeV range for sources with TS > 257 The 
red, blue, green, and black data symbols correspond, re¬ 
spectively, to 7 -ray sources detected with the Fermi-LAT 
that have been associated with FSRQs, BL Lac objects 
with and without redshifts, and blazars with optical data 
too poor to determine if the source is an FSRQ or BL 
Lac. 

From inspection of the plot, it is clear that a function 
of the form = d — k log 1^14 will provide a reasonable 
description of the data. For the entire FSRQ and BL Lac 
sample, but excluding other blazar candidates, values of 
k = 0.18 ± 0.03 and d = 2.25 ± 0.04 are deduced in the 
3LAC (jAckermann et ahl 1201511 . Comparing this value 
with the analytic expressions, eqs. ED and ED, a larger 
value of 6 is implied for SSC processes compared to EC 
processes, but in both cases consistent with 6 « 1/3. 

The typical value of 6 can also be deduced from the av¬ 
erage nonthermal blazar synchrotron SED, when fit with 
an expression of the form of Eq. (151). From X-ray analysis 
of Beppo-SAX data on Mrk 501, iMassaro et al.l (l2004f ) 
finds values of bsy ranging from 0.12 - 0.33, implying 

^ Energy flux is derived in 5 energy bands in intervals defined 
by 0.1, 0.3, 1, 3, 10 and 100 GeV. 



0.01 0.1 1 10 100 1000 


(10^'* Hz) 

pk 

Fig. 1.— Data are the > 100 MeV photon spectral index 
values r^, as a function of peak synchrotron freq uency Vs for 
blazars from the 2LAC (|Ackermann et ahl 1201 ih . Red, blue, 
green, and black symbols identify, respectively, FSRQs, BL 
Lac objects with redshifts, BL Lac objects without redshifts, 
and blazars with data too poor to determine if the source is an 
FSRQ or a BL Lac object. Left: Curves labeled by EC BLR, 
EC IR and SSC for EC processes with BLR photons, EC 
processes with IR photons and SSC processes, respectively, 
show r.y vs. Vs predictions of the log-parabola equipartition 
model using standard parameters given by eq. ED- Also, 
eo = 2 X 10“® and Uq = 10“^ erg cm“® for Lya, and eg = 
4.6 X 10“^ and ug = 10“® erg cm“® for the ~ 1000 K IR 
radiation. Thick curves give numerical calculations, and thin 
curves show analytic results, from eqs. (O and (jS)). The 
thick curves that approach constant values at large Vs are 
numerical predictions for the power-law, log-parabola model, 
eq. HH). Right: Compton-dominance Ac as a function of Vs 
for EC BLR, EC IR, and SSC processes, as labeled. The line 
with arrows has a slope of +1 in the Ac vs. Vs plane. 

a corresponding log-parabola width parameter 6 > 0.5. 
Narrow bandwidth modeling of X-ray s ynchrotron emis¬ 
sion f rom Mrk 421 gives bgy = 0.3 - 0.5 |Tramacere et ahl 
|2003), though values of bsy = 0.17 ± 0.02 (2006 15 
July pointing), bsy = 0.11 ± 0.02 (2006 April 22 point¬ 
ing), and bsy = 0.08 ± 0.03 (2006 June 23 pointing) 
are obtained in more com plete joint XRT-BAT analy¬ 
sis (jTramacere et ahl l2009ll , cons istent with a n electron 
distribution with 6 = 5bsy ~ 0.5. IChenI (j2014^ finds that 
bsy is distributed in the range 0.05 < bsy 7, 0.25, im¬ 
plying 0.25 ^ 6 < 1.25. More importantly, he finds a 
dependence of bsy on Vs, which we discuss further in Sec¬ 
tion The values of 6 deduced from spectral modeling 
tend to be larger than obtained from the slope implied 
by the spectral-index diagram. 

3.1. Standard parameters in log-parabola model 

To compare the log-parabola equipartition model with 
data, we adopt a standard parameter set, and take 

6 = 1 / 2 , 64 = L 48 = (e = Cs = Egev = 1 • (11) 

The reasoning driving the choice of the standard variabil¬ 
ity time scale is that the masses of supermassive black 
holes powering blazars—both FSRQs and BL Lacs—are 
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typically of the order ~ IO^Mq. The value ^4 = 1 or 
tyar = 3 hr corresponds to the light-crossing time across 
a size equal to the Schwarzschild radius of a ^ 10 ®Mq 
black hole, though of course shorter variability time 
scales have been recorded during spectac ular outbursts 
of B L Lac ob iec ts, including; Mrk 421 (jFossati et alJ 
[200^. Mrk 501 (I Albert et alJ [2nn7ll . and PKS 2155- 
304 (lAharonian et al.ll2007f) . not to mention the extraor¬ 
dinary VHE outburst observed with the MAGIC tele- 
sco pe from the F SRQ PKS 1222-1-216 with tyar ^ 10 
m (|Aleksic et al.l[ 20 Tl[l . The isotropic synchrotron lumi¬ 
nosity Lgyn can exceed the Eddington limit LEdd, though 
l^Edd is presumably the upper limit to the persistent ab¬ 
solute jet power (see App. B). Standard values L 48 ^ 0.1 
- 1 and L 48 ^ 10“® - 10“^ are typical of powerful FSRQs 
and BL Lac objects, respectively. At the other side of the 
time domain, tyar ^ 10 ^ - 10 ® s may be compatible with 
quiet times of blazars. 

Fig. 1 shows analytic results of Eqs. dU and dl]) for F-y 
as a function of = Vg, using the standard parameter 
set. Results of numerica l calculations, o btaine d by mod¬ 
ifying the code used in iDermer et al.l d2014al l , are also 
shown. The dimensionless photon energies for the BLR 
and IR photons used in the model are cq = 2 x 10“® 
(i.e., 10.2 eV) for BLR photons and Cq = 4.6 x 10“^ for 
warm IR torus dust emission described by an Ri 1000 K 
greybody spectrum with « 15% covering factor, giving 
an energy density of r; 10“® erg cm“®. The analytic 
results are shown by the thin lines. The numerical re¬ 
sults are shown by the thick curves. As can be seen, 
the analytic SSC and EC IR results are in reasonable 
agreement with the numerical calculations, whereas the 
analytic EC BLR results do not agree with the numer¬ 
ical results. Klein-Nishina effects already make them¬ 
selves felt strongly for target BLR photons scattered to 
1 GeV, but only weakly for target IR photons scattered 
to 1 GeV, as is clear by noting that KN effects set in at 
photon energies > TOeC®/12eo R 100 GeV for 1000 

K photons, and ^ 2 GeV for Ly a photons. The 

Thomson-regime expressions are harder than the numer¬ 
ical curves because of the Klein-Nishina softening. 

Fig. 1 also shows the effects of a low-energy power-law 
extension of the EED on the spectral-index diagrams. 
In such a power-law log-parabola (PLL P) model with a 
low-energy cutoff Lorentz factor (lYan et al.l 120131 : 
iPeuff et ai.l 1201411. the EED distribution extends eo. o 
by two parameters to take the form 

( 12 ) 

Here s is the power-law spectral index of the low-energy 
component, r is a log-parabola width parameter, and 
Vi — 'y'minl^pk- The Heaviside functions are defined such 
that H{u) = 1 when u > 0 and H{u) = 0 otherwise, and 
iL(u; a,b) = H{u — a)H{b — u). The theoretical basis for 
the form of eq. m is discussed below. Results are shown 
for s = 2 and <C 1 , in which case r ^ b, reducing the 
PLLP model to a 3-parameter model. 

3.2. Compton dominance 

The numerical results for this particular set of param¬ 
eters are seen to follow the trend of much of the data. 
Virtually no FSRQs are observed, however, with > 1. 
To obtain some insight into this, we calculate the Comp- 



(10^'* Hz) 

pk 

Fig. 2.— Same as Fig. ^ except that h = 1. Heavy and light 
downward-going curves are the numerical and analytic equiparti- 
tion model predictions, respectively, and upward going curves show 
Compton dominance for EC BLR, EC IR and SSC processes for 
the log-parabola EED, eq. Q. The thick curves approaching 
constant values at large values of Us correspond to spectral-index 
predictions of the PLLP model, eq. ns, with a —2 number index 
of the low-energy EED. 


ton dominance Ac for our model, defined here as the 
ratio of the 100 MeV - 100 GeV 7 -ray luminosity to the 
bolometric synchrotron luminosity. It is calculated from 
the relation 

, _ L..(100 MeV) 100 , 

-Tc =-—-) 

Oii/L/syn ^GeV 

where and L-y(100 MeV) are, respectively, the vL^ 
spectral index and luminosity calculated at Eq^y GeV. 
Note that a more detailed and time-intensive calculation 
would integrate the blazar SED to determine Ac- 

The Gompton dominance depends on the energy den¬ 
sity of the surrounding radiation fields. For definiteness, 
we have taken ublr = 10 “® erg cm“® and um = 10 “® 
erg cm“® in our calculations. Note that Ac scales ap¬ 
proximately linearly with uq. As Ac becomes progres¬ 
sively smaller, the corresponding blazars becomes pro¬ 
gressively less detectable as 7 -ray sources. So solutions 
should be restricted to a minimum value of Ac ■ Solutions 
should also be restricted at large values of Ac, because 
Compton drag on the jet becomes a strongly limiting 
factor, as discussed more in Section [5j Regions where 
0.1 < Ac ^ 30 may favor LSP blazars to be FSRQs, 
ISP blazars to the EC BLR, EC IR, and SSC solutions 
in Fig. [H as these values brack et measured values of the 
Compton dominance (Fig. 7 in iFin 2 [Ml). 

In Fig. (U we calculate three models in the F..^ vs. Vg 
plane corresponding to complete dominance either of Ly 
a BLR radiation (EC BLR), IR radiation from the dusty 
torus (EC IR), or internal synchrotron radiation (SSC) as 
the target photon source. Restricting the Compton dom¬ 
inance to 0.1 < Ac ^ 30 suggests that most blazars with 
jzi 4 <0.1 have 7 rays that result from scattered BLR 
radiation, while blazars with 0.1 < 1 Z 44 < 1 would have a 
mix of blazars with 7 rays made by Compton scattering 
of either BLR or IR photons, or both. At higher peak 
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Fig. 3.— Same as Figs.^and[2l except that Ecev = 0.1. 
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Fig. 4.— Same as Figs. [Jand[2l except that t 4 = 100. 


synchrotron frequencies, SSC-dominated sources would 
be most plentiful. 

The use of a different model, the PLLP EED, eq. (I12p 
with s = 2 and <C 1, is displayed in Fig. [T] and subse¬ 
quent figures by the numerically calculated spectral index 
curves that approach constant values of spectral index 
at 1/14 ^ 1. Klein-Nishina effects, described more be¬ 
low, soften the spectral index below the Thomson regime 
value of r.y = 1.5. It is interesting that essentially all 
data are softer than P-^ = 1.5, and that smaller values 
of external radiation energy density could yield typical 
measured Compton dominance values for ISP and HSP 
blazars with an external Compton y-ray component mak¬ 
ing a significant contribution to the SED. 

Figs. [5] - m show how changes in the model parame¬ 
ters affect results. Fig. [2] shows that a value of 6 = 1 


is incompatible with the combined trend of the data, 
though a values of & = 1 may be consistent with sub¬ 
populations, e.g., FSRQs. Returning to 6 = 1/2, Fig. 
[3] shows the effects of calculating the spectral index at 
EoeV = 0.1, that is, at 100 MeV rather than 1 GeV. 
Because the Compton-scattered y-ray SED becomes pro¬ 
gressively softer at larger y-ray energies, the model re¬ 
sults in Fig. [3] are uniformly harder than in Fig. [1] The 
discrepancy between the analytic and numerical results 
decreases when scattering Ly a radiation because the 
Klein-Nishina effects on the Compton cross section are 
not so great when scattering to 100 MeV as compared to 
1 GeV. The dependence on detector energy EceV should 
clearly show up in Fermi-LAT spectral index diagrams 
calculated in discrete energy ranges, e.g., 0.3 - 3 GeV and 
3-30 GeV, and should, in a statistical study, discrimi¬ 
nate between EG and SSG processes, though correlations 
between b and Vs can hide the effect. 

Fig. |4] shows how a slower variability time, with t 4 = 
100 , affects the equipartition spectral-index diagram. 
Gompared to the results in Fig. [I] the effect of longer 
variability times is to harden the spectrum. From eqs. 
([5]) and dS]), the hardening for a factor of 10 longer vari¬ 
ability time is AP^^ = -5/4 for EG processes and 
Apssc _ _35/i6 for SSG processes. Because of the dif¬ 
ficulty in measuring tyar, the variability effect on spectral 
index may be too subtle to discriminate between EG and 
SSG processes. In a statistical sample, however, more 
rapidly variable sources at equipartition would in general 
be softer, assuming that there are no underlying corre¬ 
lations between 5 and tyar, and that equipartition holds 
in the various states. 

3.3. Spectral index vs. peak Compton frequency 

F ig.[5lshows data fro m multiwavelength spectral analy¬ 
sis (jAbdo et al.ll20ldbl l of 48 bright blazars in the Fermi- 
LAT Bright AGN Sample iLBASi lAbdo et aT1l2009j) . sep¬ 
arated into FSRQs, and low, intermediate, and high 
synchrotron-peaked (LSP, ISP, and HSP, respectively, 
defined by whether iZg < 10^'* Hz, 10^'* < iZs(Hz)< 10^® 
Hz, or j/g > 10^® Hz) BL Lac objects. The upper and 
lower panels gives predictions for the dependence of P..^ 
on izc for the equipartition EG and SSG models. The 
Thomson-regime predictions, eqs. ([9]) for EC processes 
and eq. cni) for SSG processes, are plotted in black, de¬ 
pending on whether the y-ray spectral index is measured 
at 0.3 GeV (solid curves) or 3 GeV (dashed curves). The 
index is softer when the y-ray energy range used to de¬ 
termine the spectral index is larger, as noted above. 

It is worth taking a moment to explain the deviations of 
the numerical curves from the Thomson-regime expres¬ 
sions. Suppose the detector waveband EceV ^ hvc, cor¬ 
responding to the left potions of the figures for EceV = 
0.3 - 3. Consider two y-ray SEDs aligned at the same 
value of izc, one with strong Klein-Nishina effects and one 
in the Thomson regime. The SED with strong KN effects 
will be much softer at frequencies v ^ vq hy comparison 
with the one in the Thomson regime, causing the softer 
spectra when hvQ <g; EceV: sometimes dramatically so, 
compared to SEDs formed by scattering in the Thomson 
regime. 

At the other extreme hue ^ EceV, corresponding to 
the right portions of the figures, the effects of strong 
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Fig. 5. — Data points show the Fermi-LAT 7 -ray spectral index 
evaluated in the range 0.1 1- 100 GeV as a function of uFi, peak 
Compton frequency vq = blazar 7 -ray SED IIAbdo et al.l 

I2010bll . In both panels, the range to calculate the Fermi-LAT spec¬ 
tral index is at 0.3 GeV and 3 GeV for the solid and dashed curves, 
respectively, (a), upper: Equipartition Thomson-model EC predic¬ 
tions (black) are shown along with numerical predictions evaluated 
for external 1000 K radiation fields from a dusty torus (magneta 
curves) and from Lya radiation (orange), using parameters of Fig. 
[T] but with b = 1/2. (b), lower: Equipartition Thomson-model 

SSC predictions (black) are shown along with numerical SSC pre¬ 
dictions resulting from synchrotron emission with Us = 10^^ Hz 
and 10*^® Hz, as labeled. 


KN losses is to harden the low-energy portion of the 7 - 
ray SED compared to an SED formed by scattering in 
the Thomson regime (and with the same peak Comp¬ 
ton frequency). Consequently, Klein-Nishina effects will 
produce harder spectra when the detector energy range 
is less than the peak Compton frequency compared to 
Thomson scattering. 

Fig- E] shows that an EC origin in either BLR or IR 
radiation is consistent with LSP FSRQ data, but is in¬ 
consistent with an SSC origin. A similar conclusion was 
reached earlier by examining the correlation of Compton 
domina nce with core domi nance in FSRQs and BL Lac 
objects (jMever et al.ll201^ . At values of r-c ^ 10^^ Hz, 
or Ec ^ 1 GeV, the sources are all ISP and HSP BL 
Lac objects, and are compatible with either an SSC or 



Fig. 6. — Best-fit models for Mrk 501 for fixed Ce = 1 and 
letting fe vary. Data from lAbdo et ahl 1120111 1. with galactic feature 
removed. 

EC origin of the emission, given the uncertainties in tyar ■ 
In principle, however, an EC origin can be distinguished 
from an SSC origin by comparing the curvature of the 7 - 
ray component with that of the synchrotron component. 

4. NON-EQUIPARTITION MODEL FOR BL LAG OBJECTS 

Blazars may be out of equipartition, though extremely 
out-of-equipartition blazars would be less favored be¬ 
cause of the additional power required. Up to now, we 
have assumed that the equipartition parameter = 1 , 
which minimizes jet power for a given synchrotron SED 
and variability time, assuming small baryon-loading. 
Modeling of 3C 279 with Ce = 1 was possible in 
iDermer et al.1 ([2014^, though the very highest energy 
7 rays were only successfully fit by using long variability 
times with tyar ~ 10® - 10 ® s, in which case the X-ray 
emission was not well fit (cf. iHavashida et al.ll2012ll . Bet¬ 
ter hts were found in th e modeling of 3C 45 4.3 by taking 
Q between 0.6 and 3.5 (|Cerruti et al.ll20l'^ . which has a 
minor effect on the spectral slope relation.® 

It is worth asking if an equipartition situation applies 
to BL Lac objects, which would be simpler than FS¬ 
RQs by lacking significant external radiation fields. We 
apply the near-equipartition log-parabola (NELP) mod¬ 
eling technique to the 15 March 2009 - 1 Aug 2009 
multi wavelength data of the HSP BL Lac object Mrk 
501 (|Abdo et al.l 120111) . The data in Fig. [5] include 
OVRO radio observations, optical data, Swift UVOT and 
XRT data, GeV 7 -ray data from Fermi-LAT, and VHE 
data from MAGIC. Parameter values are derived using 
the Markov Chain Monte Carlo ( MCMC) te c hniqu e of 
lYan et all l|2013l) for Mrk 421 and iPeng et al.l (1201^ for 
Mrk 501, and using the 3-parameter log-parabola elec¬ 
tron spectrum, eq. Q- The fit to the TeV data is always 
bad in the Ce = 1 case. The fit with Q allowed to vary 
is obviously far better. 

The distribution of parameter values derived from the 
MCMC technique for the data of Mrk 501 is shown in 

® The fitting published in ICerruti et al.l II2013I ) lacked log- 
parabola b-dependen t factors derived in reply to the referee of 
IDermer et al.l (I2014al) . Updated values have 1 and (^s ^ 0.2. 
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Fig. 7.— Distribution of parameter values for the varying 
case. The dashed curves represent the mean likelihoods of samples 
and the solid curves are the marginalized probabilities. 



Fig. 8.— Two-dimensional probability contours of parameters. 


Fig. [3 The dashed curves are the mean likelihoods of 
samples and the solid curves are the marginalized proba¬ 
bilities.® In the fits, we run single chains and assume flat 
priors in the model paramet er spaces. Since the MCMC 
code we used in this paper (iLiu et al.ll2012t lYuan et all 
[Ml is adapted f r om C OSMOMC, we refer the reader 
to iLewis fc Bridl3 (j2002tl for a detailed explanation of 
the code about sampling options, convergence criteria, 
and statis ti cal q ua ntities. Ac c ording to th e results of 
lYan et all (l2013h . iPeng et ^ (12014 ). and iZhou et all 
(|2014l) . the MCMC method is well suited to systemat¬ 
ically investigate the high-dimensional model parameter 
spaces in fits to blazar SEDs. 

Pairs of values of (s — 2,(e — 70, and Cs = 3, Ce — 30, 
from the fitting results shown in Fig. [5] correspond to a 
change in index compared to an equipartition circum- 

® See http://cosmologist.info/cosmomc/readme.html 


stance of AF.^ = —0.2& and AF.^ = —0.146, respec¬ 
tively. Even the large deviation from equipartition causes 
a spectral-index change <0.1 unit for 6 = 0.5, and even 
less for 6 = 1/3. The typical fluid-frame magnetic field 
derived from the fits has B' 10 mC. Synchrotron self¬ 
absorption is included in the fit. Considerations about al¬ 
lowed jet power (see App. ll) restrict the departure from 
equipartition further. Thus deviations from equiparti¬ 
tion do not, on the basis of the Mrk 501 case, affect the 
spectral-index relation significantly. 

The inability of the numerical MCMC model to find 
a most favored value for tyar may reflect limitations of 
the log-parabola EED used to model the Mrk 501 spec¬ 
trum. Using a model joining a power-law at low elec¬ 
tron energies with a log-parabo la function at high elec¬ 
tron energies. iPeng et M] (120141) fit radio data down to ~ 
GHz frequencies, and obtain preferred variability times 
of tyar ~ 5 X 10^ s. 

For given values of L 48 and vn, production of the 
highest energy 7 -ray photons is assisted by going to an 
electron-dominated regime, where i/e ^ 1 and u'^ ^ u'^,. 
The larger Lorentz factor electrons required to produce 
the same value of Vs in a weaker magnetic field can Comp¬ 
ton scatter ambient photons to the highest energies. 

5. DISCUSSION AND SUMMARY 

The nonthermal synchrotron paradigm pervades think¬ 
ing in blazar physics, yet is incapable of explaining some 
of the most elementary facts, e.g., why synchrotron- 
radiating nonthermal electrons are apparently acceler¬ 
ated so inefficientl y. Rather than rea c hing values of 
«JT0r MeV ('e.g.. iGuilbert et al]ll983t Ide Jager et al.l 
Il99(if) . the peak synchrotron frequencies of FSRQs with 
i/s = 10^^ Hz are « 10^® times less than the highest 
energy synchrotron photon in the maximally efficient 
electron Fermi-acceleration scenario. Even the highest 
energy synchrotron photons from HSP BL Lac objects 
rarely exceed Ri 10 - 100 keV, orders of magnitude be¬ 
low the radiation-reaction limit. It is crucial to under¬ 
stand the reason for the low peak synchrotron frequencies 
(smaller, of course, than the maximum synchrotron fre¬ 
quency), and how they relate to source luminosities and 
SEDs, which are the basis of the blazar sequence and 
blazar divide. 

5.1. Near-equipartition, log-parabola (NELP) model 

The astrophysics developed here may point a way to 
the solutions of these puzzles by first explaining the 
spectral-index diagrams. If the radiating electrons are 
near equipartition and approximately described by a log- 
parabola EED because of the underlying acceleration and 
radiation physics, then the relationships between the 7 - 
ray spectral index F.^ and Vg and I'c precisely de¬ 
fined in the Thomson regime by functions of the form 
r.y = d — A:logi/s(c), namely eqs. dH) - (fT 0 )l . The slope 
is accurately reproduced even when Klein-Nishina ef¬ 
fects are important. Moreover, the model inputs are 
all in principle observable from near-simultaneous multi¬ 
wavelength blazar campaigns: L 48 , Vg, vq and Ac from 
spectral observations, b and Cs from SED modeling, and 
tyar from temporal analysis. As shown here for Mrk 501, 
C,e and Cs can also be deduced from SED modeling, leav¬ 
ing only the baryon-loading pu as a major uncertainty, 
which affects the jet power (Appendix iBl) . 
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The near-equipartition approach using a 3-parameter 
log-parabola EED furthermore makes quantitative pre¬ 
dictions about the dependence of observables on E-y for 
statistical quantities of blazars, or for different states of a 
single blazar. A specific example that can be performed 
with Fermi-LAT data is to determine 7 -ray spectral in¬ 
dices of a large sample of blazars of specific types, e.g., 
LSP FSRQs and HSP BL Lac objects, in adjacent energy 
bands, giving the spectral curvature. The curvature of 
the 7 -ray SED is uniquely related to the curvature of the 
synchrotron SED, depending on whether the 7 rays have 
an SSC or EC origin. The difficulty of performing this 
test, of course, is the requirement of quasi-simultaneous 
observations over a large energy range in order to provide 
a good characterization of the synchrotron SED peak and 
curvature. 

For the synchrotron spectral-index diagram, our anal¬ 
ysis shows that k = Uec = 36/4 for EC scattering 
in the Thomson regime, and k = kssc = 96/16 for 
synchrotron self-Compton (SSC) radiation. Numerical 
results show that this dependence is even valid when 
Klein-Nishina effects are important. Analysis of the com- 
bined sample of FSRQs and BL Lac objects in the 3LAC 
(|Ackermann et al.ll2015ll . k^LAC = 0.18 ± 0.03, implying 
curvatures of 6 = 0.24 if the emission arises from EC 
processes, and 6 = 0.32 if the 7 rays are SSC. However, 
it may not be correct to combine the two samples with 
different typical values of 6 in their populations. Spe¬ 
cific predictions for the slope of the P^ vs. Vs behavior, 
depending on whether the emission has an EC and SSC 
origin, should be studied for samples of blazars binned in 
ranges of 6 , because the two variables are correlated, with 
larger curvatures, 6 ft! 1 for FS RQs, compared to 6 < 0.5 
for BL Lac objects (|Chenll20T3l . Insofar as the SSC com¬ 
ponent seems less dominant in FSRQs (Cs ~ 0.2) than in 
BL Lac objects (^s ~ 1), the effect of this correlation on 
the spectral-index diagrams also has to be considered. 

In principle, underlying correlations of P-y with tyar 
can be examined with the increasing number of simulta¬ 
neous multiwavelength blazar SEDs. Limitations of the 
log-parabola function to describe the EED remains a cen¬ 
tral assumption that can be relaxed, though not without 
associated theoretical or numerical efforts. 

5.2. Departures from equipartition 

One of the uncertain parameters is the electron 
equipartition parameter Q, here defined as the ratio of 
nonthermal electron and positron energy to magnetic- 
field energy throughout the volume of the radiating re¬ 
gion. Assuming — l,Cs = 1 gives the model predic¬ 
tions shown in the Figs.[T]~|4l As shown in Adp.IbI large 
departures from equipartition are not allowed if the ab¬ 
solute jet power is required to be less than the accretion 
power, which in turn is assumed to be bounded by the 
Eddington luminosity.^® From the results of App.jB] one 
possibility is that the SSC bolometric luminosity in the 
SEDs of large Compton-dominance FSRQ flaring events 
should be small compared to the bolometric synchrotron 
luminosity (that is, fs ^ 1 ) for compatibility with sub- 
Eddington jet powers. 

Spectral modeling of the FSRQs 3C 279 and 3C 454.3 

Counter-examp les to this assumption are claimed 
(IGhisellini et al]|2014l') . See Add. IbI 


is poss ible for Ce — 1 (jCerruti et al.ll20f^ iDermer et al.l 
For the BL Lac Mrk 501, a large departure from 
equipartition is required to get a good spectral fit, as we 
have shown, but even in this case, the effect from this out- 
of-equipartition condition on P.^ is small. The deviations 
from equipartition giving the best fits to the SEDs of Mrk 
501 show that large Ce, electron-particle-dominated fits 
(with correspondingly weak magnetic fields) are favored 
to fit HSP BL Lacs extending into the TeV regime. 

5.3. Extensions of the log-parabola model 

The log-parabola function, eq. o, is motivated by 
second-order Fermi acceleration theory where MHD tur¬ 
bulence in the emitting fluid systematically accelerates 
particles to form a curving EED (see Section 2). An 
equally compelling scenario combining first- and second- 
order processes considers a power-law distribution of par¬ 
ticles injected downstream of a shock into a turbulent re¬ 
gion where second-order processes broaden the distribu¬ 
tion, so that the EED approximates the PLLP function, 
eq. (Il 2 |). 

The full PLLP model has 5 parameters, but we have 
treated in Figs. 1 - 4 the important case of an EED with 
a —2 number index extending to low energies without 
cutoff. This EED makes a low-energy boundary to the 
spectral-index data near the Thomson value of P-y = 3/2. 
Remarkably, this is as hard as the hardest Fermi-LAT 
blazar spectral indices measured. So if Compton scat¬ 
tering is responsible for the formation of the 7 -ray SEDs 
of HSP blazars, as is undoubtedly true for the bulk of 
the radiation, then the PLLP model would give a simple 
explanation for the lack of blazars harder than P-y = 3/2. 

The situation is complicated, however, by the possi¬ 
bility that if the EEDs had low-energy cutoffs rather 
than power-law extensions to low energies, then HSP 
blazars in the LAT band would tend to be dim and hard 
to detect. So the apparent lack of blazars harder than 
r.y = 3/2 could be a selection effect rather than a limit 
imposed by the radiation physics. Searches in the Fermi- 
LAT 7 -ray data for blazars harder than P.^ = 3/2 would 
test whether a PLLP model is preferred over a LP model; 
the detection of such hard blazars would rule out the form 
of the PLLP considered here. 

Except when i/g 10^® Hz, the existence of a low- 
energy power law in the EED makes only a small differ¬ 
ence to the SSC predictions compared to the pure LP 
model. Figs. 1-4 show that the SSC predictions tend 
to be slightly softer than the data. A low-energy cutoff 
in the EED could harden the SEDs, making it possible 
to attribute an SSC origin to the 7 -ray components to 
all HSP blazars. The discovery of blazars harder than 
r.Y = 3/2 in the Fermi-LAT energy range would support 
this interpretation. 

5.4. Blazar types in the P.^ vs. plane 

We now ask why there are essentially no FSRQ blazars 
with i/g ^ 10^^ Hz. The answer is likely to involve 
the dynamics of increasingly higher synchrotron-peaked 
near-equipartition jets which, when finding themselves 
in an external radiation field, are subject to a radia¬ 
tion force opposite to the direction of motion that acts 
on the nonthermal electron population. In the ideal 
one-zone model considered here, there is no radiative 
drag from synchrotron and SSC processes, only from 
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EC processes (|Trainacere et al.]l2011f) . The larger val¬ 
ues of (5 d and 7 ^^, for increasing Vs implies a correspond¬ 
ingly larger radiative drag when external radiation fields 
are present that would eithe r slow the jet plasma down 
(|Ghisellini fc Tavecchidl^OlOH or prevent it from reaching 
such large T factors in the first place. 

Rather than treating the jet dynamics, which is beyond 
the scope of the present enquiry, we quote simple analytic 
expressions for the synchrotron and external Compton 
SEDs from log-parabola distribution, eqs. m and 0 , 
which effectively answers the question of how the Comp¬ 
ton dominance (oc radiative drag) grows with increasing 
Vs- From eq. 0 , to keep the apparent synchrotron lu¬ 
minosity constant requires ^ constant. Eq. Q 

shows that the apparent external Compton component 
grows oc 6^uo"fpi.. The ratio of the EC to synchrotron 
component is essentially the Compton dominance, which 
goes ~ 5‘^uq/B'‘^ oc vi 4 ^l]J^uq. So the Compton dom¬ 
inance and radiation drag grow oc j^i 4 , all other things 
being equal. This confirms the behavior shown in Figs. 

which deviates at large values of z^i 4 due to Klein- 
Nishina effects. Depending on baryon loading, the jet 
could be quenched before escaping the BLR, making an 
unusual, short flaring event. In such an unstable situa¬ 
tion, persistent emissions with large z /14 in dense external 
radiation environments might not be possible. 


5.5. Maximum particle energy 

The near-equipartition log parabola model can be used 
to derive expressions for maximu m escaping pr oton or ion 
energy Emax- Starting with the iHillaa (jl984fl condition 
in the form Emax = ZecB'S^tyar implies 


Emax(eV) = 1.4 X 


, ^1/4 pl/8 

/ ^4 xi/8 Jl J 2 _ 


(14) 

using the dependences given in Table [T] for the NELP 
model. For a BL Lac object with L 48 < 0.01, the only 
way to accelerate ultra-high energy cosmic-ray protons to 
> 10^° eV occurs when <C 1, that is, in a magnetically- 
dominated jet. It is interesting to compare this expres¬ 
sion with the formula 


Emax ( 6 ^) 


= 2 X 10^°Z 


VeB(Lp/t/ 10 ^® erg s ^)/ee 

f/lO 


(IWaxmanI l2004t iFarrar fc GruzinovI l2009t 
iDermer fc Razzaoua |20101 l which was Sso derived 
from the Hillas condition, where Lph is the isotropic 
bolometric photon luminosity, and cb and Ce are the 
fractions of jet power going into magnetic field and 
electrons, respectively. 

An interesting feature of a combined lepto-hadronic 
blazar model using log-parabola functions for the particle 
distributions is that the synchrotron radiation-reaction 
limit for protons is ~ 200r GeV, a factor mp/me greater 
than the electron limit. The evolution of the combined 
lepton synchrotron/SSC and proton synchrotron SEDs 
with would favor a proton synchrotron component in 
the same large magnetization (Ce 1 ) regime where the 
electrons are incapable of making high-energy radiation. 


5.6. Blazar sequenee and hlazar divide 


Multiwavelength data from any given blazar display 
a rich array of spectral and variability properties. The 
spectral properties of blazars in this analysis are reduced 
to i/g, VQ, r.^ and 6 , while the variability properties 

are reduced to tyar- 

The spectral-index diagrams show robust correlations, 
which we explain as a consequence of relativistic blazar 
jets with different powers and in different environments, 
within which are entrained relativistic electrons that can 
be described by log-parabola EEDs. By relating the syn¬ 
chrotron peak frequency and synchrotron SED, which 
mirrors the EED, to the spectral index of the 7 -ray SED 
formed through EG or SSC processes, the dependence of 
T-y = d—k log Vs is easily derived in the Thomson regime. 
Moreover, the equipartition relations imply specific pre¬ 
dictions for underlying correlations. 

The spectral-index diagrams are one side of a trian¬ 
gle relating T.^, Vs (or Vc), and The term can 

either be the apparent isotropic synchrotron, 7 -ray, or 
total bolometric luminosity. The other two sides of the 
triangle are vs. Vs or vc, the blazar-sequence rela¬ 
tions, and r.y vs. Lj, the blazar-divide relation. 

Our work illuminates one side of the triangle, namely 
r.y vs. Vs or vc- Regarding the blazar divide, suppose 
as a first approximation that the typical mass of a su- 
permassive black hole is 10 ®Mq, then Fermi-LAT data 
shows a significant change of spectral index at the Fermi 
divide of Ly = 10"^® erg s“^. If the apparent 7 -ray lumi¬ 
nosity is 10 % of the apparent jet power, and the beaming 
correction is ^ 100 , th en the divide is at L/Lbaa = 0.01 
l|Ghiselhni et al.l[2009ll . and this would also represent the 
Eddington ratio below which the external radiation field 
energy density becomes small. 

Extremely weak dependences, if any, are seen in the 
r.y vs. Ly blazar-divide plots on either side of the di- 
vide. Within b l azar subpopulations (see Fig. 39 in 
lAckermann et all 1201 111 , the dependences of T.^ on Ly 
are also weak. Other than near the divide itself, there 
is no clear dependence of r.y on blazar luminosity. In¬ 
deed, any such dependence is predicted to be weak, as 
can be seen from eqs. 0 and 0 , which show that T.^ oc 
—blogLsyn/8 for EC processes, and T.^ oc 61ogLsj^„/32 
for SSC processes. 

To explain the blazar sequence relating and Vg or 
vc requires jet physics outside the scope of the present in¬ 
vestigation. Rather than saying why blazars of a certain 
type can exist, however, we can suggest why blazars dom¬ 
inated by EC emission require low-synchrotron peaks. 
The presence of any appreciable external radiation field 
would produce a Compton drag that decelerates the bulk 
flow or prevents such a near-equipartition situation that 
would produce a synchrotron SED peaking at such large 
Vs from forming. 


6. SUMMARY 

To conclude, we have used an equip a rtition blazar 
model ing approach (jCerruti et al.l l2013t IDermer et al.l 
l2014a|l to explain the correlations of Fermi-LAT 7 -ray 
number spectral index r.y with peak synchrotron fre¬ 
quency Vs and peak Compton frequency vc- This ap¬ 
proach assumes a one-zone model fit to the broad¬ 
band emis sion, so that emissions from, e.g., e xtended 
VHE jets (jBottcher et al.ll2008t lYaii et al.l 1201211 . spine- 
sheath structures ( Ghisellini et al.l l2005lL decelerat- 
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ing jets (|Georganopoulos fc Kazanai l2003lj . or VHE 
emissions induced by UHECRs produ ced by the jet 
(|Essev et al.l 120101 : iTakami et al.l 120131) . are assumed 
not to affect the 7 -ray spectral indices or peak frequen¬ 
cies. Within this framework, the trends in the spectral- 
index diagrams are reproduced in a model with equiparti- 
tion conditions and a log-parabola electron distribution 
with b = 1/2. This conclusion holds even for out-of- 
equipartition conditions limited by absolute jet power to 
be sub-Eddington. 

The broadly distributed data in the spectral-index di¬ 
agrams suggest that a better model comparison would 
consider a distribution of parameter values to define a 
preferred model region in the T.^ vs. Vg and r'c diagrams. 
Such an approach de pends on kno wing whether the cor¬ 
relation of b with i^s (iChenl 1^0141) is robust, if one is to 
sample from a distribution in b values. Nevertheless, al¬ 
lowed regions in the spectral-index diagrams in Figs. 1 - 
5 are already defined by the heavy solid curves, depend¬ 
ing on whether internal SSC or external EC BLR or EC 
IR processes dominate the formation of the 7 -ray SED. 
A distinct trend in the boundaries on the spectral index 
diagrams in Figs. 1 and 3 are found for the PLLP model 
that can be tested with Fermi-LAT analyses in different 
energy ranges. 

The r.y vs. Vs boundaries defined by the dominance of 
internal SSC or external EC IR or EC BLR processes 
have, furthermore, a very different shape when 6=1 
(Fig. 2) compared to 6 = 0.5 (Fig. 1). This can be tested 
by subdividing the Fermi-LAT 7 -ray spectral indices in 
different ranges of 6 . Whether the log-parabola function 
or the PLLP model, eq. m with a low-energy electron 


index s = 2, better approximates the EEDs can be tested 
by searching for Fermi-LAT sources with P.^ < 1.5. 

The weak dependences of T-,, on changes in Lgyn found 
in eqs. ([5]) and ([ 8 ]) are consistent with the weak de¬ 
pendences of 7 -ray spectral index T.^ on on either 
side of the blazar divide. A physical explanation for 
the change of the radiation environment of blazars at 
Ri O.OlLEdd, though a reasonable model assumption, 
would make sense of the blazar divide. 

This leaves open the blazar sequence relations, which 
can ultimately only be understood from the physics oc¬ 
curring in the magnetospheres of the supermassive black 
holes powering the blazars. Near-equipartition blazar 
synchrotron sources with Vg ^ 10^^ Hz would suffer in¬ 
creasingly strong radiation pressure in an environment 
with dense external radiation fields, which could explain 
the absence of HSP FSRQs. Supermassive black-hole 
jets are most luminous when their emissions are coolest, 
that is, when their peak synchrotron and Compton fre¬ 
quencies are lowest. The near-equipartition log parabola 
blazar model provides a constrained system that explains 
the spectral-index diagrams, and points to studies that 
could allow for a deeper understanding of the blazar se¬ 
quence and blazar divide. 


The work of C.D.D. and J.D.F. is supported by the 
Chief of Naval Research. We thank Dr. Matteo Cerruti 
for discussions about spectral fitting, and the anony¬ 
mous referee for constructive questions and the recom¬ 
mendation to consider the PLLP model. 


APPENDIX 

A. (5-FUNCTION THOMSON-RECIME SSC DERIVATION WITH LOG-PARABOLA EED 


We derive the form of the vL^, SED for the SSC component in the Thomson regime assuming a log-parabola func¬ 
tion of the EED and employing (5-function approximations for synchrotron and Thomson scattering. The comoving 
Thomson-scattered synchrotron self-Compton spectrum for isotropic distributions of photons and nonthermal rela¬ 
tivistic electrons is given by 




^3J2 


d'y' 



m ') Kii) 


da{e) 

de[ 


(Al) 


Here, yi' is the cosine of the angle between the directions of the interacting electron and photon, e = 7 e'(l — fj!) 
is the invariant collision energy, and dai^jdd-^ is the differential scattering cross section. We use the (5-function 
Thomson scatterine cross section (i(Tfe)/(ie' = — v'^e'll —u')l fea. iG.ddlilDermer & Menonll2009D . From ea. o, 

— K'y ^ /ipk^ photon spectral density = e'L'(e')/ 47 r/oi?(,^e'^meC^, where /o is a geometry 

factor, = cSntvar, and tyar{^ + z) is the measured variability time (|Dermer et al.ll20f4^ . For the synchrotron target 
photon spectrum, e'L'^y„{e') = €Lsyn{e)/Sfj = /Sf^. Plugging these expressions into eq. (lAll) . and using the 

5-function to solve the fj! integral, we find 




utvK' 


'^D47r/oi?f7pfeepfc Ji/r, 


dy y 


-6-blogy 




dx 


-6-blogx 


(A2) 


Here xg = y/Ajy, where A = e'i/ 27 p|,e/j. = e/cpk^ssc- The interior integral can be solved by noting, to good 
approximation, the logarithmic term is slowly varying compared to the x~^ term. The value of this integral is then 
y^-5-b\ogxe After some manipulations, we obtain 




2aTvK'j'yi^ i_|iog^ 

bnfoR'^ V 26 


(A3) 


Comparing with eqs. m and o shows that the SSC spectral index is given by eq. @ with 6 replaced by 6/2 and epk 
by epk,SSC, leading to eq. dH). 
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B. JET POWER IN THE NEAR-EQUIPARTITION LOG-PARABOLA MODEL 
We consider jet power with the addition of baryons and photons. The baryon-loading factor 77 ^; = and 

where is the fluid energy density of both thermal and nonthermal protons and ions, and u'^ is the 
nonthermal lepton energy density, including both electrons and positrons, For convenience, the the rmal electron and 
positron energy density is assumed small. The absolute jet power for a two-sided jet is given by (iCelotti fc FabianI 
II 993 I : ICelotti fc Ghiselli'nill2008t iGhisellini fc Tavecchidl201^ 

Ljet = 27rr^^/3r^cu^/[1 -|- <^e(l + ^w)] + Lph , (Bl) 

where the absolute photon power Lph comprises synchrotron and SSC radiations, each assumed to be emitted isotropi¬ 
cally in the jet frame, and EC radiations, with its comparatively narrower b eaming (l]Permei]ll995D . The absolute photon 
powe rs depend on the observing angle 9 through the Doppler factor Jp (jChisellini fc Tavecchid 120101 : iDermer et al.l 
[ 2 OII), giving the absolute jet power in the NELP model: 


Ljetieig s-^) = (1+iV^)" {4.0x10 


-.44 

\ 



+ \/G{^ + Vbl)] + 


2{L 


ISO 

SSC 


T ISO \ 

' ^syn) 




o r 4SO 

/ 

5^d 


(l+N^n (B2) 


where the observing angle 9 = Nr/T 1 and P » 1. In this expression, L^ssc^ ^ec measured 

apparent isotropic bolometric synchrotron, SSC, and EC luminosities, respectively. The other two terms in eq. (IB2I) 
correspond to the magnetic-field, oc l/y^, and the particle power, oc ^^e(l -|- iju)- The additional factor of (1 -l-iVp)^ 
narrows the focus of the 7 -ray beam. Eq. (IB2I) can be rewritten as 

.5/8 A/2 r - ^ 

Ljetierg s"^) = 4.0 x 10^^ (1 + \ J(i + (-+«;), (B3) 

\JOSs) ^ ^14 ^ ^ 


where 


w = 


Ceji+ r]bl) 
1 T Vph 


(B4) 


and the radiation loading 

r,ph=^= 5.4/oG [1 + Cs + 0.6Aec(1 + N^)^] . (B5) 

Ub, 

The isotropic bolometric photon luminosity + L'^ec + -^SSC = where the external Compton 

dominance Aec = ^EcZ-^syn- 

From eq. (IB3I) . the minimum power condition is defined by the condition w = 1. When the baryon loading factor 
r]bi ^ 1 and the radiation loading rjph <C 1 , the minimum power condition is defined by = 1 - If the baryon-loading 
is arbitrary, but rjph <C 1, the minimum power condition is defined by (e = 1/(1 + When rju ^ 1, the minimum 
power condition also corresponds to a highly magnetized jet (in terms of the electron energy density), with a larger 
jet luminosity by a factor y/1 + r]bi at minimum jet power compared to a pure electron/positron jet. When rjbi and 
r]ph take arbitrary values, the minimum power condition is defined by Ce = (1 + Vph)/i^ + Vbi), and the minimum jet 
power incr eases oc r]bi){l -I- ijph). 

Eq. (IB3I) gives the absolute minimum power to make the observed radiations from a blazar jet. For example, if 
the angular extent 9j of the jet exceeds l/L, the power is increased by ~ (LS^)^. When observing at 0 ~ l/L, the 

minimum jet power, Ljet ~ 3 x lO'^^L^g® erg s“^ can only be increased by one to to orders of magnitude before 
exceeding LEdd ~ 1-3 x lO'^’^Mg erg s“^ for a 10 ®Mq black hole, unless one demands unusually high photon efficiencies. 
For HSP BL Lac objects, with Lspn ^ 10"^® erg s“^, there is no great difficulty in satisfying the Eddington limit, even 
far from equipartition. However, these very same objects are believed to be accreting at a rate ;< O.OILedDi so even 
in this case, large departures from equipartit i on ca nnot be tolerated. 

Based on Fermi-LAT data. IChisellini et al.l (|2014f) argue that the absolute jet power Pjet is larger than the accretion- 
disk luminosity, which is approximated as 10 times the BLR luminosity. They also approximate Pjet ~ lOProd, with the 
absolute radiation power P.ad ~ where the factor kf = 8/3 for synchrotron/SSC processes and kf = 32/5 

for EC processes, and the bulk Lorentz factor L of the radiating jet’s plasma outflow is stated to be in the range 
10 < r < 15. The un derly ing assumption is that the observer is looking at 0o = l/L to the jet axis. Inspection of the 
photon power in eq. (IB3I) shows how uncertain this assumption is given how much brighter fluxes are along the jet 
axis compared to fluxes from sources at 9q « l/L. Using the relation 1 -|- Ap = 2r/(5D, the photon power in eq. (IB2|) 
for a 2 -sided jet is 

= (1 + «?)- + (1 + Kf (B6) 

When viewing down the jet axis, these powers are ft! 16 (syn chrotron/SSC) and r; 64 (for EC) times less than the 
values used in the expression for Prad bv IGhisellini et ahl (120141) . The most powerful 7 -ray sources have the largest core 
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dominances and brightn ess temperatures ie.g.. iPushkarev et alJl2009t iKovalev et alJl2009HLi et al.ll20l4 however, see 
Savolainen et al.l (l20ir)li i. suggesting that these sources are also the ones viewed almost along the jet axis. Indeed, 
Jorstad et al.l (|2005ll . Fig. 25, find no blazar with viewing angle 0q > 2/r, whereas a large number of both BL Lacs 
and FSRQs have 0o < l/2r. A severe overestimation of the radiation and therefore jet power is made by not taking 
this effect into account. 

For very powerful FSRQs like 3C 454.3, which has L 48 ~ 1, a curious feature arises. Great flares exceeding 
^ with large Compton dominance > 100 can be allowed while maintaining absolute jet power 

Ljet ^ Afidd only if 0 <C 1/F and Cs 1, that is, implying a small SSC component relative to the 

synchrotron component. The effect of decreasing Cs is to increase Jd and narrow the Doppler cone, making the beaming 
factor even smaller, so that extreme apparent EC 7 -ray powers lead to absolute jet powers that are sub-Eddington. 
Spectral modeling to give the relative SSC and EC powers depends on X-ray observations, for example. Swift and 
NuSTAR, Fermi-LAT observations at GeV energies, and ground-based VHE air Cherenkov arrays. 
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